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Abstract 

In this paper we investigate a system describing electrically charged 
particles in the whole space R^. Our main goal is to describe large time 
behavior of solutions which start their evolution from initial data of small 
size. This is achieved using radially symmetric self-similar solutions. 
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1 Introduction 

We analyze the Debye-Nernst-Planck system of the form 
ut = + V ■ (ijV0„), 

2 (1^) 

(t)u = —E2 * II in R , 

where u(x,t) : x R+ R and -£2(2) ~ ^ log \z\ is the fundamental 
solution of the Laplace equation in R^. We supplement this system of 
equations with the initial condition 

u{x,0) = uo{x) in R^. (lb) 

One may think about this system as of a single parabolic equation with 
a nonlocal nonlinearity. 

The history of studies of this model is quite long. The system was 
introduced by Nernst and Planck in the nineteenth century, and then 
modified by Debye in twenties of the twentieth century (see 02] )• It de- 
scribes the evolution of the density u{x, t) of electrically charged particles 
interacting through the potential (f)u- We encounter systems of this type, 
also for several unknown densities, in the electrochemistry, theory of semi- 
conductors and physics of plasma [U [H [16] [20] . 
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Many difFerent methods have been used for the analysis of this system 
since that time (see [H UJ SI [S] and references therein). However some 
problems still remain unsolved. 

Similar problems where the first equation of (|lap is replaced by (jju ~ 
E2 * u, appear in the chemotaxis theory and they are used for describing 
the evolution of particles which interact via gravitational forces [51 1141 
115) . Some reasonings from that theory can be adopted without essential 
changes but questions related to the long time asymptotics have com- 
pletely different answers. For example, Banach type theorems giving the 
local existence can be applied in both cases (see the next section). On the 
other hand, we know that gravitational problem does not have solutions 
for large initial masses (see [24] and the references therein), which is not 
the case of electric forces. 

Whenever it would not cause any confusion we omit time or spatial 
variable. Moreover (/ * g) [x) denotes the usual convolution, i.e. (/ * 
g) (x) — Jjj2 f{x — y)g(y) Ay, and S{t) is the heat semigroup given by 

{S{t)u) (x) = f G{x- y, t)u{y) dy, (2) 

where G{x,t) = e~'^' is the Gauss kernel. Here and subsequently 
||m||p denotes the usual L''(R^) norm of the function u{x,t), with respect 
to the spatial variable. C is a generic constant which can vary from line 
to line. The Schwartz class of rapidly decreasing functions is denoted by 
and its dual by ^'(R^). 5o{x) denotes Dirac delta function at the 

point 0. 

In each dimension if u{x,t) and 4>{x,t) solve the system (|la|l then 
)<^u{\x, \^t) and (f>{\x, \^t) also do. As we will see the total charge of 
particles, namely the integral J^2 u{x,t)dx, is preserved in time. How- 
ever, solutions invariant under the above scaling preserve the charge only 
in K'^. So that we expect essentially nonlinear asymptotics of generic solu- 
tions. Moreover, self-similar solutions which describe long time behavior 
decay faster than the Gauss kernel when \x\ 00. The two dimensional 
case differs much from the higher dimensional case, where the long time 
behavior is described simply by the heat kernel [4]. However, in |4] au- 
thors use a completely different approach involving energy and entropy 
dissipation methods. 

First we prove (following methods used in [3| [17]) the existence of 
solutions for the system (|lap . Then we study self-similar solutions in 
radially symmetric and general setting. The last section is devoted to the 
large time behavior of solutions. 
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2 Existence of solutions for singular ini- 
tial data 



We are mainly interested in mild solutions framework. We look for solu- 
tions of the integral equation 

t 

u{t) = S{t)m + j S{t-T)V ■ {u\7(j)) (t) dr, (3) 



where 4'{x,t) and uoix) are as in ([l}. Two major problems connected 
with this approach appear. One is: in what sense the solutions of the 
integral equation (|3]) satisfy the differential problem (|la|l . It is obvious 
that the classical solutions of (|la|l satisfy the equality (|3]). Unfortunately, 
the reverse implication is not always true (see also a discussion of this 
problem in [XT]). 

The second problem is: in what sense the initial condition is fulfilled. The 
answer for this question will come with the definition of the space in which 
we solve the integral problem (|3}- 

Let us notice here that the term V- {uV(l>u) can be written as B{u, u) = 
V ■ (uVcpu) where B(u, v) is a bilinear form defined by 

B{u,v) = V ■ [uV(j>v), with — — -B2 *v, 

for u, V from some Banach space X, which in our case is a space L'^^^{M^^) 

(see a discussion before a proof of Lemma [2.11 below) . 

We introduce here the notation of the rescaled version of a given function 

/(^) 

where A is some real and positive parameter. In the case of the function 
depending on space and time, this kind of scaling (and, therefore, the 
notation) is used only with respect to the spatial variable. 

By a simple calculation we show that the bilinear form B(u, v) has 
scaling order equal to (cf. 17 ) which means that for every A > and 
every u,v £ X 

B{ux,vx) = {Biu,v))^. 
We use frequently very well known estimates of the heat semigroup 

\\s{t)f\\^ct'/--'/^\\f\\, , 

(4) 

i|V5(t)/||^^CT^/--^/'-^/=||/||, , 

where 1 ^ q ^ p ^ 00, C — C{p, q) is some positive constant, t > and 
f{x) £ L^(R^). Let us also recall the weak Young inequality 

\\g*h\\-,^C{(3,-f,g)\\h\\^ . (5) 

for all h e L"(R2), \g{x)\ <: \x\-'^l'' and i/q + i//3 = 1 + V^- For the proof 
see, e.g., [191 Remark (2) in Theorem 4.3]. 
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The inclusions ^(R^) C L^/^{R^) C J^'(R^) are continuous and the 
norm || • II4/3 is translation invariant. Moreover the estimation 

\\S{t)B{u,v}\l^^ ^ Ct-'^'\\u\U^,\\v\U^, (6) 

holds, for every u,v £ I/*/^(E^). 

Indeed, using the second inequality in (|4]) with p — '^/s and q = 1 and the 
Holder inequality, we get 

We use then inequality ([5]) with 7 = 4, /3 = 2 and a — ^jz to obtain 

||u||4/3||V(^„||4 = ||lt||4/3|||a;r^ * ^ ||lt||4/3|b||4/3, 

which concludes the proof of (|6} . 

It means that the space I/*''^(R^) is adequate to the problem ([l]) in the 
sense of [171 Definition 3.1]. We construct a homogeneous Besov space 
modeled on the Banach space X — I/*''^(R^), as 

BA" = |u G ^'(R^) : supf'''''||5(t)u||^^^ < c5o| , 

and we introduce a norm on this space by 

||u||ba- = supf'''''||5(t)wjL, . 
t>o ' 

■ 1/2 

The space BX is in fact the classical homogeneous Besov space B^j^ ^. 

■ 1/2 

The norm || • ^bx is equivalent to the original norm in Blj^ ^ introduced 
by the dyadic decomposition. For a further discussion on this topic we 
refer the reader to [171 Example 4.2] and the references therein. 

As it was in [171 Lemma 5.2] let us prove here the a priori estimate 
for the bilinear form B{u, v). 

Lemma 2.1 There exists a constant C > independent of u, v and t 
such that 

\\Smi.U,v)\\^^ <: Ct-^/^|lw||4/3 ||«||4/3 , (7) 

for all u,v £ X . 

Proof We begin with the fact, that there exists some function C{t), 
independent of u and v such that 

=sups^/*||5(s)5WB(n,t;)|| s;C(t)|h||4/3||H|4/3 

for every t > Q. The proof consists of two steps. First let us assume that 
< s < 1. Then also < s^'^ < 1 and because ||5(t)/||4/3 II/II4/3, the 
inequality (|6]) yields 

s'^^S{s)Sit)Biu,v)\l^,^ \\Sit)Biu,v)\l^^ < Ct-'/^||w||4/3lb||4/3 ■ 
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e 4Ai! 



In the case s > 1 we proceed in a similar way but first we use the semigroup 
property S{s)S{i) = S{a -\-t) = S{t)S{a) to obtain 

s''^\\S{s)S{t)B{u,v)\l^^ = s'/'\\Sit)S{s)B{u,v)\l^^ 

^ S^^''\\S{S)B{U,V)\\^^^ ^ CS"'/"" ||?i||4/3 ||u||4/3 ^ (^11^114/3 ||w||4/3 • 

Now, we show that C{t) is some positive constant. First let us check the 
scaling property 

S{t)B{ux,vx) = {S{X^t)B{u,v))^ , (8) 

for all u, V e L^/'-'iR^) and for all real A > 0. For the proof let ^ G .y'(R^). 
Then 

(9) 

As we mentioned before the scaling order of the form B{u, v) equals zero 
which yields 

S{t)B(ux,vx) = S{t){B{u,v)) ^ = (s{\h)B{u,v)) ^. 
Then we use scaling properties 

IIMI4/3 = A-^''^||/||4/3 and WfWsx = X^WhWsx , 
to show that 

\\S{\h)B{u,v)\\^^ = \^\\{S{\h)B{u,v))^\\^^ = X^\\S{t)B{ux,vx)\\^^ 
^ \^C{t)\\ux\U/:, ||«a||4/3 = C(t)\-^\\u\U/.,\\v\U/:, . 

Now we fix f = to > and we put A = ('/to)^''^. Then we obtain the 
estimate 

\\s(t)B(u,v)\\^^ ^ c{u)tfr'/^\\u\U/,\\v\U/,, 

which concludes the proof. ■ 

The space where we are looking for the solutions of the integral equa- 
tion ([3]) is defined as 

^' = ^([0,00) -.BX) n 

n {f{t) : (0,00) ^ L'/'iM?) : suvf/'\\f(t)\U < 00}, 

(10) 

and we equip this space with the norm 

ll/lla: =max|sup||/(f)||fl;,, supf^/*||/(t)||4/3}. 

t>0 t>0 J 



The space X with the above norm is a Banach space. We understand 
^([0,oo) : BX) as the space of -BAf-valued measurable functions which 
belong to ([0, oa) : BX) and attain the initial condition of the problem 
((T| as t — > 0"'" in the sense of tempered distributions, namely 

J <t>(x)u(x,t) dx J (j)(x)uo(x) dx, 

for all (l> € as t goes to 0+. 

Now let us formulate and prove the main result of this section. In the 
proof we combine the reasonings presented in [3] and |17l 118) . 

Theorem 2.2 There exists an e > such that for all uq G BX satisfying 
\\uo\\bx < £, there is a global in time solution u{x,t) of the problem ([T]) in 
the space pO|) . This solution is unique among all the functions u £ BX 
satisfying the condition sup^^Q t''^\\u(t)\\x < 2e. 

Proof We show that the operator 

t 

J^{u){t)=S{t)uo + j S(t-r)B[u,u)(T)dT 



is a contraction in the set 



B. 



•,^{ue.T: Mt)\\.^ s; 2e] 



for some small £ > 0, to be specified later. 

In the following estimates we use frequently the elementary inequalities 
such as the Holder inequality, the heat semigroup estimates Q and the 
estimates of the bilinear form B{u, v) from © and (0 



^ l|uo||g_y +C ^supr^/''||M(r)||^^^^ J {I - wy'^^^w'^^^ dw 



and similarly 

t'%yriu){t)-^{v)it)l^^ 
t 

=^ Ct'/' Jit- r)-'^%u - «)(r)j|,/3 (lhMlL/3 + \Hr%,) dr 



^ Cesvipt'^^\{u-v){t)\\ ■ 
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In the same way we deal with the estimate 



2 

dr 







and 

\\.yr{u)it)-,yr{v){t)\\^^^ Cesupt'^'Wiu^vmll^^ . 

Combining the above inequalities and choosing e > such that Ce^ < e 
and Ce < 1, we obtain 

IKWII,*; < 2e , 

\\{u-v){t)\\^ , 

for every u, v £ I/''''^(R^). It follows that ^ : is a. contraction 

in Bg. Thanks to the Banach fixed point theorem we obtain the existence 
of a mild solution of the problem ([1]). This solution is unique in the set 
Be and is global in time. ■ 

Since we are particularly interested in self-similar solutions, it is de- 
sired to know what does the condition ||uo||ba" < £ mean for the initial 
condition being a distribution. Namely, let uq{x) = MSxq{x). Then 

4/3 



/ 2 \ 

E2 ^ K2 / 



e 4t \ ^ / g 3 4t 

ATvt J (47rt)V3 J ATvt 

M^/^ 1 f 1 , 

e ' dz. 



3 (47rt)V3 tV4 y 47rt 



R2 



Then 



limsupf'''*||M5(t)5^J|4/3 ~ M*''', 

which means that atomic measures are admissible as initial data (|lb|) but 
only with small masses of single atoms. 

Moreover let us notice, that if uo{x) is a function from L^(R^) - it 
does not contain singular part - then 

limsupi'"^*||uo||4/3 — 0. 

Above shows that the part lim sup^^g t^''* || • II4/3 of the norm || • || a: measures 
only the singular part. 
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3 Local and global existence in L^(IR^) space 



Using the same framework of mild solutions we prove the local in time 
existence of solutions of the problem ((3} in the space C([0, T]; L^(R )) n 
C{[0,T]; L^CR^)) for some p > 1 and an arbitrary initial condition uq G 
Z/^(R^) nI/''(R^). Moreover, we show that these solutions can be extended 
to global in time solutions. This situation differs from the case where par- 
ticles interact through gravitational forces and through the chemotactic 
attraction. The solutions in that case blow up in L^(R^) spaces for suffi- 
ciently large initial data (see [141 115) ). 

In the following we use the same inequalities as before, i.e. the heat 
semigroup estimates (|4]) and the weak Young inequality ([5]). 

Let us define the space 

Xt = C([0,r];L'(R^)) nC([0,r];L''(R^)), 
which is Banach space with the norm 

\\u\\xj, = sup ||u||i -I- sup \\u\\p . 

Proposition 3.1 Let ^ ^ p < i. For any uo £ L^(R^) n L''(R^) there 
exists some T > and a unique solution of the integral equation 

t 

u{t) = S{t)uo + j S{t- r) V ■ (uVfli) (r) rfr, 



such that u{x,t) £ Xt- 
Proof Let 

t 

.yr{u)(t)=S{t)uo + j VS{t-T)- {uV(j}){T)dT. (11) 


Then Q implies that 

t 

lp(t) = ||._/r(n)(t)-<S(t)uo||p j {t^TY'--''--''^\\{uV(p){T)\\^AT, 



for every 1 ^ p, q ^ oo. We use the Holder inequality to obtain 

t 

Tp{t)^C f (t-r)^/--^/'-^/^||u(r)|| ||(V0)(r)||^dr. 

J ^ p-q 



We can exploit ((Sjl with a = pq^2p-2q 'f^ ^ ^ and 7 = which yields 

t 

Ip{t) f (t-r)^/'-V--^/^||u(r)|| ||K(r)|| dr. 

J " pq-2p-2q 
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We would like to have pq^2p-2q ~ P which is possible for q = and all 

4/3 ^ p < 4. Consequently, we have 

sup \\^{u){t) -S{t)uo\\p 



0<t<T 



<C( sup Ht)\l^ j{t-Tr'^-dT 



2 ^ 

^CT'-''-( sup \\u{t)\\ \ [{l-sy'/'ds, 





where the last integral is finite since —3/4 < — i/p ^ —1/4. 
Similarly, we obtain 



sup^||^(w)(t)-5(t)wo||i 



2 1 

CT^-V^f sup \\uit)\\ ] [il-s)'-'/'ds, 



with the same constraint for p, namely ^3 ^ p < 4, which yields the 
inequality —1/2^1 — 2/^^1/2. 
An analogous calculation gives 



supJjK{u){t)-J'{v){t)\\p 

sup ( sup ||w(t)|| + sup \\v{t)\\ 



^ C-yi-Vp 



0^ 

and 



sup \\J^{u){t)-J^{v){t)\U 

^CT^-'/" sup f sup \\u{t)\\ + sup 

for all w, w € At and Ys ^ p < 4. Without loss of generality we can 
assume that < T < 1. Then 

|K(«)WL, ^ \\u,\U,+CT'-'''-\\u{t)\\^^ 

and 

\\^{u){t)-^{u)[t)\\^^ 

^ cT^-''-\\u{t) - v{t)\\^^ (ihwlL^ + Ht)h,) ■ 
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For arbitrary _R > we can choose T > so small that the operator ,yY 
is a contraction in the ball Br{0) C Xt- The Banach fixed point theorem 
gives the existence of a local in time solution in the space Xt- From the 
construction this solution is unique in the ball Br{0). But the uniqueness 
is a local property and the above solution is also unique in the space Xf 
for r > possibly smaller than T. Indeed, let us suppose that there 
is another solution in the space Xt contained in some ball Bf^{0) with 
R < R. Taking T < T we construct the unique solution on the time 
interval [0,T] but in the ball _B^(0), which contradicts the uniqueness of 
the previous solution. ■ 

Proposition 3.2 Let 4/3 ^ p < 4 and uo{x) e L^(K^) n L''(R^), // 
additionally uo{x) is nonnegative, then for allO <t <T, u{x,t) ^ and 



luix,t)dx = luoix)dx. 



R2 

Proof The conservation of positivity for the solutions of the problem 
((TJ is a consequence of the analogous property for the system on bounded 
domains. Indeed, there has been proved in [2j Proposition 1 and Propo- 
sition 2] that solutions with initial conditions ^ uo{x) £ L^{^1) n L^{fl) 
are nonnegative for t > 0. The proof involved Stampacchia truncation 
method. The same approach permits to prove nonnegativity of the solu- 
tion u{x,t) either by approximating u{x,t) on the sets {|a:;| < A^} where 

G N, by solutions u{x,t) = Uk{x,t) of the problem posed on the sets 
< A'^ -|- A:} where A: £ N or by a direct rewriting of the Stampacchia 
scheme (see also |12l Lemma 1] and Theorem I, Lemma 4.1]). 

The L^(R^) norm of the solution is properly controlled, as it was jus- 
tified in the proof of Proposition 13.11 namely 



sup \\Af{u){t)\\i^\\uo\\i+CT^-'^''\\u{t)\\^ . 

0^t<T 

Integrating the Duhamel formula (|lip over we get 

t 

^ u{x,t)Ax = J S{t)uo{x)dx + J J \/S{t-T)- {u\/(j}){T)dTdx. 

But the last integral on the right hand side is equal to zero, since we can 
apply the Fubini theorem, and the first one is equal to Jj^j uo{x) dx. ■ 

We obtain the global in time existence of the solutions using the fol- 
lowing a priori estimate. 

Lemma 3.3 Let u{x, t) be a solution of the problem ((!]) with nonnegative 
initial condition uq{x) £ L^{M?). Let M = Jjj2 uo{x) dx. For every p G 
(1, co) such that 2'^~^ < p < 2^ for some fc G N, the inequality 

\\u{;t)\\v ^ Cik)Mt-^^-'M 

holds with a constant C{k) given explicitly. Moreover 

\\u\\oo ^ CMt-^ 

holds for some constant C > 0. 
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We base the proof of the above lemma on the following interpolation 
inequality 

Lemma 3.4 Let 1 ^ q,r ^ oo, j, m G N, ^ j < m and a £ [^/m, 1]. 
Moreover 



- = ^ + a +(l-a)-, 

p n \r n y q 

where n is the space dimension. Then 

J2 \\D-u\\,^cf J2 WD-uWrY \\u\\l-\ (12) 

|a|— |a|— m 

This inequality was introduced by Gagliardo-Nirenberg-Sobolev and some- 
times is referred to as the Nash inequality 21 . 
Proof of Lemma [33] 

Step 1. Let p — 2. Multiplying the first equation in (|lap by u and 
integrating by parts, we get 



^ ~ J l'^"* ~ uAu — uV ■ (uV<^)) da; 



1 d 

(13) 



^^Il«ll2 + IIV^II^ + ^Jvu^- V<^da; 



= i^lhlli + l|Vu||^ + - / u'dx. 

I 

From Proposition 13 . 2 1 the solution u{x,t) in for all t > 0. Then we 
obtain the inequality 

^lhlli + 2||V«||i«:o. (14) 



Then taking a — ^2, q — m — 1, n = r — p = 2, j = Oin (|12p we get 

\\uh^C\\Vu\\T\\u\\T (15) 



which is equivalent to 



u\\t ^ \\Vu\\l (16) 



CM, 

where Mo = ||u(0)||i. Using and {TG]) we get 
Integrating the above inequality yields 

II /'j-M|2 ^ 1 ^ CMq 1 



which ends this part of the proof. 

Step 2. p = 2^ . We prove it by recurrence. Let as assume that 



rlfc— 1 rlfc — 1 /i-,fc— 1 



sC afc_iMo" t-'^' (17) 
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where 



] '~^h for k — \, , . 

' '^2'- 2(a^_j)2 for fc>2. ^ ^ 



In the inequahty p5|l we set instead of u which yields 

(||n||^:)'<C||Vu^'"||^(||u||^ti)'. (19) 
The same calculation as in (|13|l with the multiplication by ~^ gives 

^rf+^^^l|V«^^"llU0. (20) 
Combining together inequalities p7|) . p9|) and (|20|l we get 

l|2''-\2 



C(||n|| 



2'=-l - 



d „ ..o'' 2*= - 1 



|. II2* , Z -1 V||M||2fc, 



\ 2 



< ^ r 



{\\u\\llf ^ CA'/f 2'=-2(a,_i)^ 

Let us denote f{t) = \\u(t)\\^t. and integrate the above inequality over the 
interval {0,t). We obtain 

1^1^ 1 



fit) /(O) " CM2'"2*-2(afc_i)' 



1 1 2'«-l 1 



/(t)^:C2'=-^(afe_iMo^')V(2''-^) 

which gives 

WuWll ^akMtt-^^''-^\ (21) 

Step 3. General case. For each p £ (l,oo) there exists some k £ 'N 
such that 2''"^ < p < 2*. Then it is sufhcient to use inequality (I2ip and 
very well known interpolation inequality 



with a= Then 



1-1/2* 



At the end let us notice that 



(l/p-l/2'')/(2'=-l) , ^l/p-l/2fe , 1/2 
k ■ k k 



Step 4. p = oo. In the last step we use the property 
\\U\\oo = limsup \\u\\p. 

p — > oo 
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Then it is enough to show that lim^^^oo ra^f^ < oo. But at — (7"'= 2™*=, 
where 

-1 for fc = 1, 

2wk-i + {k-2) for fc > 2. 

and Wfc = 2* — 1. It is not difficult to see that for all k £ N, Wk < 0. Then 

Um aV^ — lim — — —r- < C < oo 

which ends the proof . ■ 

As a simple corollary from Proposition 13. II and Lemma 13.31 we obtain 
the main result of this section, namely 

Theorem 3.5 Let 4/3 p < 4 and uo{x) £ ^^(R^) n L^{M?). Then, 
there exists the unique solution of the integral problem ([3]) in the space 
C([0, oo) : L^{9?)) n C([0, oo) : L^{R^)) . 

4 Radially symmetric self-similar solutions 

As we mentioned in the introduction the system (|lap is invariant under 
the scaling 

\^u{\x,\^t), 4){\x,\^t), (22) 

where < A G R. 

In this section we assume that u(x,t), <j){x,t) are self-similar, i.e. they 
are invariant under the scaling ((22} . The most important from our point of 
view are the solutions beginning their evolution from the initial conditions 
of the form M5o{x). Such an initial condition is obviously radial. Thus we 
assume from now on the radial symmetry of the solutions u{x,t), (j){x,t). 
Following the integrated density method (see [31 [3 [25]) we introduce the 
new variable 

r 

Qir,t)^ I uix,t)dx^2.lsuis,t)ds, 

where Br{0) is the closed ball of radius r centered at the origin (0,0). 
Integrating the system (|lap over such balls and performing simple com- 
putations, we rewrite the problem (|la|l - (|lbp in the form 

Qt — Qrr Qr "T QQr^ 

r Iivr 
Q(0,t) = 0, for t^O 

lim Q{r,t) = / uo{x)dx = M, for t ^ 0. 



Using self-similarity we can again change the variables: Q{r,t) = 2n^{y) 
for y = r'^/t. This way we transform our problem to the following one 

e"(y) + i?'(j/)-^C(y)C'(y) = 0, (23a) 
m = 0, ?'(0)=a. (23b) 
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Formally, the boundary conditions implied by (llbfl are 
C(0) = and lim ^{y) = M. 

But it is more convenient to investigate equation (|23a|l with the boundary 
conditions H23b|l for some positive real a. As we will see, both problems 
are equivalent. The reader can find similar results considering the problem 
(|23p (existence of solutions, regularity, etc.) in 7 but methods used there 
are different. 

The main result we prove in this paper is the following 

Theorem 4.1 For all M > 0, there exists a unique solution of (|23ap . 
satisfying the initial conditions ^(0) =0 and limy^oo ^(y) = M . 

First we justify some a priori properties of the solutions of (|23|l . Then 
we are able to construct a suitable space where the Banach fixed point 
theorem can be applied. 

Lemma 4.2 Let ^{y) be a solution of (|23|) with y G [0, Yq)- Then 

i) C'(0) > 1/2 implies i'{y) > 1/2, for all y G [0, Yq); 

li) ifO< C'(0) = a < 1/2 then < ^'(y) < a and -i/s < ^"(y) < for 
all y G [0,yb]- Moreover, in such a case < ^{y) < ay on the whole 
existence interval; 

Hi) if (,{y) is global and < ^'(0) < 1/2, then limy_oo exists; 

Proof (i) In order to compute ^"(0), we let y tend to 0^ in the 
equation (I23al) . Then 

r(o)^c'(o)Qm-i)^f(«-i) 

sines 

C'(0)= hm iM^. li^ iM. 

H-.0+ y y-,o+ y 

In the case a > 1/2 we have £,"{y) > for some right neighbourhood of 
y — 0. Now let us define 

yo = sup [y G [0, Yo) : ^'{s) > \ for all s G [0, y)}, (24) 

and suppose that we have ^'(yo) = V2 ^ point yo < Yq. 

It follows immediately that ^'(y) is decreasing on some subinterval 
of [0,yo). Since for y G [0,yo) we have ^(y) > iy, then by (|23ap . the 
second derivative (y) is positive on the interval [0,yo). It follows that 
^'(y) cannot decrease on this interval, which contradicts our assumption. 
Thus, we have C'(y) > V2 for all y G [0, Vo)- 

(ii) It is easy to check, that if < a < 1/2 then ^'(y) < 1/2 on [0, Yq). 
It is enough to note that ^"(y) < 0, which implies that the first derivative 
^'(y) is strictly decreasing, so ^'(y) < ^'(0) = a. It remains to prove that 
^'(y) is positive. Integrating ^'(y) < '^/2, we get that ^(y) < y/2 for y G 
[0, Yo)- Assume on the contrary, that there exists a point y G [0, Yo), such 
that 5'(y) = 0. Then ^"(y) = 0. We have to consider two possibilities: 
^- i' {y) < Oj ■S. > 0; both, for y in a small, right neighbourhood of 

y, i.e. [y, y + e) —: TZ^iy), with sufficiently small e > 0. 
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ad 1. We know that ^'(j/) is decreasing in TZt{y). On the other hand, since 

m^f(,)(f*-i)>o. 

and < y/2 we deduce that (,"{y) > for y > y. Then the 

function ^'(y) cannot decrease for any y > y; 

ad 2. Here £,'{y) is increasing in 7it{y). But in the same manner as in the 
previous case, we see that ^"{y) < 0. 

A contradiction ends this part of the proof since integrating the above 
inequahty and using the initial conditions we obtain the remaining part 
of the conclusion. 

(in) We have shown that S,'{y) > and £,"{y) < 0. We have to 
consider two possibilities: either the limit of £,'{y) at 00 is equal to or it 
is a positive constant less than 1/2. We will exclude the second case. 

First of all we rewrite (I23al) as 



e'jy) _ fgy) 1 

Using (ii) we show that the right hand side of the above equality is strictly 
negative, namely 

^^^^S = -Ul^a) <0. (25) 



2V2 

We can write (I25p in the form 

(inC'(y))' ^ -S, 
and after the integration from to some y > yo, we obtain 
InC'(y) <lnC'(0)-5j/, 
C'(y)< e'"''e-*^ 
0<C'(y) <ae-'^ 

Let J/ go to cx). Then 

lim ^'{y) ^ a lim e"''^ ^ 0. 

y — ►oo y — >oo 

Integrating again the inequality £,'{y) < ae"*" on the interval [0,y], we 
obtain 



y 



~Ss 1 

e as, 



e(j/)-?(0) <a / 

(26) 



5(y)< 7(1-6 ") = - — (l-e''^ ' 



which concludes the argument. ■ 

To prove the existence of solutions of the problem (|23|) we transform it 
into an integral one and apply a Banach fixed point theorem. Multiplying 
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equation (|23a|) by e"''*, integrating and using the initial conditions (|lb|) . 
we obtain 



"2^ 



^(y) = 4a(l-e-''/^)+i | (e"''"/ ^C(^)C'(s)d.) dt. 



We look for fixed points of the operator 



H(C)(j/)=4a( 



1 - e 



-t/4 



— C(s)C'(s)ds)dt, 



in the space 

T = C^[0,Fo]n{C: C(0) = 0, C'(0)=a, < ^'(y) < i, for j/ G [0, Fo] } , 
endowed with the usual norm 

iieiir= sup sup i'{y)^m\c + m\c^ . 

Lemma 4.3 There exists lo > such that for a given < a < 1/2 the 
problem (|23[) has a solution in the space T. Moreover this solution is 
unique. 

Proof We widely use here the properties of ^{y), which follow from 
the choice of the space T and Lemma [4.21 Let us begin with 



||H(C)||c<4a( 



1 - e 



^ 4a 1 - e 



^ 4a 1 - e 



-Y-o/4 



^ 4a(l- e-''«/^ 



+ 1 sup / e-'/^ / ^\^(s)^'(s)\dsdt 



+ i sup / e-'/* / ^|^(,,)|||e|Ieidsdt 



V t 

+ lM\\h sup /e-'/^ / e/'dsdt 



O^y^Ya 
u 

+ 2||C||ci sup {y + 4e 



-y/4 



+ 2yo||ci 
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Similarly, we have 



||H(C)||ci = sup 



ae +--( / fe 









a+i sup e-"/^ / ^|e(.)||C'(s)|ds 



^'i+ollCIIci sup e ^'^ 1 as ds = 



^a + 2a||C||ci sup (l _ e""/^) = a + 2a(l - e"''"/*) ||C||ci 
which gives 

\\H{i)\\T ^ a + 2a(l - e-"«/*) (2 + ||C||ci ) + 2yo||Cllci • 
To show the contraction property for the operator Ti we estimate 

s/i 



\\niO-n{v)\\c<^ sup f e-'/' f ^\as)^'{s)-v{sW{s)\dsdt 



y t 

< i sup / e-'/^ / — \as)\ ■ IC'(s) - V(s)|d.dt 



+ i sup f e-'^' f ^\as)~v{s)\-Wis)\dsdt 



y t 

^ — sup / / ae'''* sup |^'(s) — 7;'(s)| ds dt 

2 o^j/ssro J J o^s^yn 



/ i e=/* sup |e'(s)-r;'(s)|dsdt 



1 

+ 9 sup 



||e-r;||ci(2a + l) sup fy + 4 e-"/* - 4 

< y;,(2a + l)||C-»7llci • 

We estimate the term ||7i(^) — 7i(?7)||ci in the same way, to get 

\\m)-nmr^C{Y,)\\C-^\\r. 

Choosing a sufficiently small Yq we justify, as in Theorem 12.21 the exis- 
tence of the solution to the integral problem. To complete the proof we 
have to show that ^(y) constructed above solves the differential equation 
(|23ap . For all y £ (0, Yo) our solution is a smooth function since the sin- 
gularity in the coefficient Yh disapears. Then we can consider the same 
equation but on some interval not containing zero. The uniqueness of the 
solutions ends the argument. ■ 
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Lemma 4.4 The solution obtained in Lemma \4-3\ exists for all real y G 
[0,oo). 

Proof Let us assume that the solution ^(j/) is not a global one. Then, 
there exists a finite 

Vmax ~ sup{y : ^(y) exists and is in the space T} . (27) 

Using monotonicity and boundedness of the solution on the existence in- 
terval given by Lemma 14.21 we can let y tend to Kmax- Let us denote 

max 

) by a and ^ (ymax) by /3. Then 

e"(?/) + ie'(y)-^C(j/)C'(y) = o, (28) 

e(ymax)=Q, C'(ymax)=/3. (29) 

is an ordinary differential Cauchy problem which is not a singular one on 
the interval [Ymax, Yg) for some Yq > Kmax- Then there exists the solution 
S,{y) of this problem, on the interval [ymax,^^) (see, e.g. [13]). Since the 
function 

'ay) for ye[o,Yo), 
J{y) for yelYoX) 

belongs to T on the interval [0, Yq), we get a contradiction with the max- 
imality of Ymax- B 



v{y) = 



Now we can prove the main result of this paper. 
Proof of Theorem 14.11 Lemmas 14.31 and 14.41 assure the existence 
of the unique solution of the problem (I23ap - (|23b|l for y £ [0, co) and 
a € [0,1/2). Moreover, from the inequality (I26|l we obtain the estimate 
from above for total mass of the solution, i.e. ,^(+00). 

As in [231 Theorem 21], we use the comparison principle for a second 
order ordinary differential equation, to show the estimates from below of 
our solution. In order to do this we have to regularize our problem and 
show that all the estimates are independent of the introduced regulariza- 
tion parameter. Indeed, let solve the equation 

F{a) := aJ(y) + ^'M - Uy) ay) = o, (30) 

with the initial conditions ^e(O) = 0, Ce(0) ~ 9 for some positive e and 9. 
Then for h{y) = A{l - e""^), with A = 

F{h{y)) = - ^,'Ae-^y ~ ^^e-™ (^^^(l ^ ^~'') - \) 



fiA 



'My 



4(2y + e) 



2y + s 
(2y + e)(4M-l)+4^(l 



which is not positive for all y £ R^, uniformly with respect to £. Existence 
of the solutions to the problem (|3Up is a standard fact (see, e.g., [13]). It 
suffices to show that family {^£}o<e<i is equicontinuous on every finite 
interval [0, Yo]- Then, passing to the limit e ~> we conclude the proof. 
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Now we see explicitly the equivalence of the both formulations of our 
problem: with the initial conditions ^(0) — 0, C'(0) = o, and the boundary 
conditions ^(0) — 0, (,{oo) — M. Taking ^'(0) = a we solve our system 
with total mass -r^^ and if we want to know solutions with mass M we 
have to take ^'(0) = in (f23b)l . 

5 Asymptotics 

Before proving the main result of this section we recall one more technical 
fact (for the proof see [171 Lemma 6.1]. 

Lemma 5.1 Let h(x) be a nonnegative function in L^(0, 1) such that 
Jgh{x)dx < 1. Moreover, we assume that a{t), fi{t) : (0, oo) R"'" 
are hounded and such that 

1 

a{t) ^ fi{t) + j h{x)a[xt)dx, (31) 



for all t G (0, oo). Then limt^oo Pit) — implies limt^oo ci(t) = 0. 

In the case of the problem ((T} we are able to prove analogue of [181 
Theorem 3], namely the following theorem. 

Theorem 5.2 Under the assumptions of Theorem \2.Sl let u, v be two 
solutions of the problem ^ with the initial conditions uq, vo £ BX, 
respectively, obtained in Theorem \2.^ Then, for sufficiently small e > 
the relations 

lim t^'*\\S{t){uo~vo)\\^,=Q (32) 

and 

lim t'^^\\u{-,t)-v{-,t)\\^,^^Q (33) 

are equivalent. 

Proof The integral representation of the solution ((3| yields 



i{t)-v{t) = S{t){uQ~VQ) + j S{t - t){B{u,u) - B{v,v)){t) At 



t 

= S{t)(uo ~ vo) + j S{t~ t){B(u,u~ v) - B{v,u~ v)){t)At. 



Then 



+ ^/(*-^)"'^'(lhM|l4/3 + II^Mll4/3)l|(^-^)Mll4/3dr. (34) 



According to assumptions in Theorem 12.21 



supf"'''''j|u(t)|| ^ 2e and supt^^''||i;(t)|| , ^ 2e. 
t>o t>o 
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Multiplying the inequality (|34|l by t^^"' and denoting 

a{t) = t'^'\\u{t) - v{t)\l^^, m = t'/'\\sm^o - ^o)||,/3, 

we get 

t 

a(t) s; /3(f) + c J{t- tY''^H^i\~^i^^ 



X (rV^||n(r) 11^/3 + rV^||«(r) 11,^3) II (n-^,)(r) 11^/3 dr 
1 

= /3(f) + 4eC / (1 - u;)"^''''w"'/^Q(it;f)du.. 



Since = 4£C(1 — £ 1/^(0,1) and, if necessary, at the 

expense of decreasing of £ > 0, the quantity 4eC /g^(l — 'u})~^^'^w~^^^ Aw 
is strictly less then 1, we can apply Lemma [5. II to obtain that the relation 
^ implies 

To obtain the reverse implication we proceed in a similar way. First let 
us estimate 



1 



(35) 



Then we notice that 



sup {wty^'^\\u{wt) -v{wt)\\ ^ sup t^^'^\\u{t) -v{t)\\ 
™e[o,i] ' t>o ' 



^supi' \\u{t)\\ +supt' \\v{t)\\ <oo, 
t>o ' t>o ' 



and use the Lebesgue dominated convergence theorem. We obtain that 
the right hand side of the inequality psp converges to zero as t tends to 
00. Then we deduce that the relation (|33|l implies p2p . ■ 



We derive as a corollary of Theorem 15.21 a remark on local asymp- 
totic stability of self-similar solutions or, in other words, the self-similar 
asymptotics of solutions to (|3]) with small initial data. 

Corollary 5.3 Suppose that uo satisfies the assumptions of Theorem \2.Sl 
Wo (a;) dx = M and vo{x) = Um{x,1), where Um is the unique self- 
similar solution with Um{x, 1) dx = M . Then 

lim t'^''\\u{-,t)-UMi-,t)\L,^^0 

holds. 

Proof Indeed, if 110(2;) is a regular initial condition with ^^0(2;) dx = 
M, then 

lim^t^^'^\\S{t){uQ - vq)\\^^,^ =0, 
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cf. |17] . because the L''-asymptotics of the heat semigroup is determined 
by the Gauss kernel multiplied by M, namely 

Um t^''^''\\S{t)vQ - A/5(t)5oj|p = 0, 1 p oo. 

This is the reason why we took the self-similar profile at the time t = 1, 
not the singular data at t = 0. ■ 

We expect that the local result in the above corollary can be extended 
to the global asymptotic stability of self-similar solutions on all L'' -spaces, 
that is 

lim t^-'^''\\u{-,t) ~UMi-,t)\\ =0 (36) 

t— >oo P 

(which is in fact implied by (|36p for p = 1), for each uq G BX and Um 
with M — uo{x) dx. 

The proof of this result seems to be beyond the scope of methods in that 
paper. Indeed, imitating the variational approach as was in [S] for the 
problem with the gravitational interactions (in particular: energy and 
entropy functionals as well as entropy dissipation relation), one may first 
prove a priori estimates on u and then a bound on u{x, t) log t/^^^'i) dx, 
leading to (f36)) . 
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